Quantum Annealing with Jarzynski Equality 



Masayuki Ohzeki'', Hidetoshi Nishimori 

"Department of Systems Science, Kyoto University, Yoshida-Honmachi, Sakyo-ku, Kyoto 606-8501. Japan 
''Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro-ku, Tokyo 152-8551, Japan 



Abstract 

We show a practical application of an well-known nonequilibrium relation, the Jarzynski equahty, in quantum com- 
putation. Its implementation may open a way to solve combinatorial optimization problems, minimization of a real 
single- valued function, cost function, with many arguments. It has been disclosed that the ordinary quantum com- 
putational algorithm to solve a kind of hard optimization problems, has a bottleneck that its computational time is 
restricted to be extremely slow without relevant errors. However, by our novel strategy shown in the present study, we 
might overcome such a difficulty. 
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1. Introduction: Quantum Annealing 

To reduce power loss in electric circuits, we have to 
minimize the circuit length. This kind of problems are 
formulated into a more generic task to minimize or max- 
imize a real single-valued function of multivariables, 
cost function. This is called optimization problem |[l|. 
To solve these problems is one of the most important 
tasks and has broad applications in science and engi- 
neering. Well-known examples with discrete variables 
are satisfiability problems, exact cover, maximum cut, 
Hamilton graph, and traveling salesman problem. 

Most of the above exemplified cases are mapped 
into a generic problem to find the ground state for a 
kind of systems seen in statistical physics, spin glasses. 
Its Hamiltonian is denoted as //q in the present study. 
One of the generic algorithms to solve optimization 
problems in reasonable time by exploiting resources in 
physics is quantum annealing (QA) iSSHl- In QA, we 
introduce artificial degrees of freedom of quantum na- 
ture, noncommutative operators in order to induce quan- 
tum fluctuations. 

H(t)^f(t)Ho + {l-f(t)}Hi, (1) 

where /(f) is assumed to be a monotonically increasing 
function satisfying /(O) - and /(r) - 1 . The anneal- 
ing time is denoted by r. The transeverse-field operator 
for spin glasses is often used as a quantum fluctuation. 
Hi - - Yji'^r The quantum annealing starts from the 
ground state of //] , a uniform linear combination of all 
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spin configurations on the basis of cr~.. The quantum 
system is driven by gradually decreasing the quantum 
fluctuation according to /(f). The adiabatic theorem 
guarantees that we can reach a nontrivial ground state 
of Hq after sufficiently slow quantum sweep r^. ~ 1 /A^, 
where A is the energy gap of the instantaneous quantum 
system (HJ L5]- If we consider the cases in which the 
quantum system as in Eq. ([1]) has a minimum energy 
gap vanishing as A ~ exp(-Q'A^) for increasing the sys- 
tem size A^, QA does not work weU in reasonable time 

m. 

2. Jarzynski Equality 

To overcome the above difficulty in hard optimization 
problems, we propose a novel method in conjunction 
with a theoretical piece from non-equilibrium statisti- 
cal physics, the Jarzynski equality (JE), in the present 
study |l8i|9t]. The Jarzynski equality is written by an 
well-known expression as, 

where the angular brackets denote the average over all 
realizations in a predetermined process starting from an 
initial equilibrium state and W is the work done during 
the process. The partition functions for the initial and fi- 
nal Hamiltonians are written as Zo(y8) and Zrifi) with in- 
verse temperature fi, respectively. We here shortly recall 
the formulation of JE for the classical system on a heat 
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bath. Let us consider a thermal nonequiUbrium process 
in a finite-time schedule fo = < f < f„ = nSt. Ther- 
mal fluctuations can be simulated by the master equa- 
tion. The probability that the system is in a state cr^. at 
time tk is denoted as PCcrj.; f^). The transition probabil- 
ity per unit time 6t is defined as M(cr|(.+i|o-A.; f^). In the 
original formulation of JE, the work is defined as the 
energy difference merely attributed to the change of the 
Hamiltonian, but we can construct JE also in the case of 
changing the inverse temperature by defining the work 
as -/3W{tk) = -{J3{tk+x) - l3{tk))Ho for the state cr. The 
left-hand side of JE can then be expressed as 

(e-^»^) 

n-\ 

^ ^ j~j |g-/5H'(/l)g<5fM(cri^,|crs;/,,)| 

xP(cro; fo), (3) 

where ^(ctq; fo) denotes the initial equilibrium distribu- 
tion. Even if the transition term exp(5fM(cr^.+i|cri.; ?^)) 
is removed in this equation, JE is trivially satisfied as 
one can simply confirm. A non-trivial feature of JE is 
in the insertion of the transition term. From Eq. (O, 
it is straightforward to prove JE by use of the detailed- 
balance condition. An observant reader may think the 
above formulation without any consideration of quan- 
tum nature is not available for the application to QA. 
Nevertheless we can apply the classical JE to QA by aid 
of the classical-quantum mapping ifioll . 

3. Classical-quantum mapping 

The classical-quantum mapping leads us to a special 
quantum system, in which the (instantaneous) equilib- 
rium state of the above stochastic dynamics can be ex- 
pressed as the ground state. A general form of such a 
special quantum Hamiltonian is given as Hq{cr'\a-, t) — 
I - e'^"«M(cr'\cr; t)e.^^"<''. This Hamiltonian has the 
ground state as I^PeqCO) = y/Z{ij. The 

ground state energy is 0, which can be explicitly 
shown by the detailed-balance condition. On the other 
hand, the excited states have positive-definite eigenval- 
ues, which can be confirmed by the application of the 
Perron-Frobenius theorem. 

In the above special quantum system, we can deal 
with a quasi-equilibrium stochastic process as an adi- 
abatic quantum-mechanical dynamics in QA. Let us 
consider QA for the above special quantum system by 
setting the parameter corresponding to the temperature 
r — > oo (j6 — > 0). This condition gives the trivial ground 
state for Hi in the preceding section with an uniform 



hnear combination, similarly to the ordinary QA. If we 
change T — > very slowly, one can obtain the ground 
state of Hq, which expresses the very low-temperature 
equilibrium state for Hq, the cost function of the opti- 
mization problem that we wish solve. We however con- 
sider to construct a protocol with the same spirit as JE 
by using the special quantum system to overcome the 
bottleneck of the ordinary QA as proposed in the fol- 
lowing section. 

4. Quantum Jarzynski annealing and its application 

We prepare a trivial ground state with a uni- 
form linear combination as the initial condition in 
the ordinary QA. This initial state corresponds to 
the high-temperature equilibrium state I'I'eqCfo)) 
exp(-y8(ro)//o/2)|o-) with j6(fo) <«; 1. We introduce flie 
exponentiated work operator Wexp(f/(:) = exp(-(j6(ft+i) - 
/3(tic))Ho/2). It looks like a non-unitary operator, 
but we can construct this operation by considering 
an extended quantum system |lll [l2|] . If we apply 
Wexp(tk) to the preceding quantum state \^eq(tk)}, it is 
changed into a state corresponding to the equilibrium 
distribution with the inverse temperature fi{tk+\). Af- 
ter then, the time-evolution operator U(o-'\o-,tk+\) - 
exp(-i6tHq{o-'\o-, tk+\)/ h) also does not alter this state, 
since it is the ground state of Hq{o-'\o-, tk+\). The result- 
ing state after the repetition of the above procedure is 

n-l 

mtn)} n {Wexp(f*)t/*+i(crt+ik,; f,)} . 

Xl^eq(fo)) (4) 

This is essentially of the same form as Eq. ([3]). We 
measure the obtained state by the projection onto a 
specified state cr'. The probability is then given by 
Kcr'|*I'(f„))p, which means that the ground state we wish 
to find is obtained with the probability proportional to 
exp(-y6(f„)//()). If we continue the above procedure up 
to /3(t„) » 1, we can efficiently obtain the ground state 
of Hq. This is called the quantum Jarzynski anneal- 
ing (QJA) in the present study. Most of the readers, 
who are familiar with the ordinary computation, have 
considered that it may seem unnecessary to apply the 
time-evolution operator U ((Tk+i |cr^; ft), which expresses 
change between states by quantum fluctuations, at the 
middle step between the operations of the exponenti- 
ated work operators Wexpih)- The time-evolution op- 
erator does not mean an artificial control but describes 
the change by quantum nature during quantum com- 
putation, which is inherent property in quantum com- 
putation. However we here remember the nontrivial 
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state in a short annealing and implies that we may over- 
come the difficulties in hard optimization problems and 
solve them in a reasonable time. The present result is 
nothing but preliminary one. We should address the 
problem on practical efficiency for several interesting 
hard problems we wish to solve in the future study LI3I1 . 
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Figure 1: The performance of QA (triangles), and QJA (cir- 
cles). The reference curve describes the instantaneous Gibbs- 
Boltzmann factor. 



point of JE. Even if quantum nature affects the instan- 
taneous quantum state, the property of JE guarantees 
that we keep the quantum state expressing the instan- 
taneous equilibrium state. If one considers to simu- 
late this procedure in "classical" computers, we have 
to need the repetition of the pre-determined process to 
deal with all fluctuations in the nonequilibrium-process 
average, since we assume an ensemble in equilibrium 
in the formulation of JE. However, when we implement 
QJA in "quantum" computation, we operate QJA to a 
single quantum system in principal since the classical 
ensemble is mapped to the quantum wave function. We 
do not need the repetition of the same procedure differ- 
ently from the classical case 11211 . 

Let us take a simple instance to search the minimum 
from a one-dimensional random potential, which is for- 
mulated as the Hamiltonian Hq - -Yj'iLi Here 
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Vi denotes the potential energy at site / and chosen ran- 
domly. We employ a linear schedule for tuning the pa- 
rameter /3 from to 100. Figure [T] shows the plot for 
QJA (upper circles), which are fixed along the refer- 
ence curves (solid curve) representing the instantaneous 
Gibbs-Boltzmann factor. In contrast, QA (lower tri- 
angles) can not sufficiently find the ground state since 
we consider a very short annealing is considered in this 
case. 
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5. Summary 

We consider an application of JE to quantum com- 
putation as QA to solve the optimization problems by 
using the classical-quantum mapping. As we expected, 
this protocol keeps the quantum system to express the 
equilibrium state for the instantaneous inverse temper- 
ature. The result by QJA shown here gives the ground 



3 



